Spin squeezing in optical lattice clocks via lattice-based QND measurements 
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Quantum projection noise will soon limit the best achievable precision of optical atomic clocks 
based on lattice-confined neutral atoms. Squeezing the collective atomic pseudo-spin via measure- 
ment of the clock state populations during Ramsey interrogation suppresses the projection noise. 
We show here that the lattice laser field can be used to perform ideal quantum non-demolition 
measurements without clock shifts or decoherence and explore the feasibility of such an approach 
in theory with the lattice field confined in a ring-resonator. Detection of the motional sideband due 
to the atomic vibration in the lattice wells can yield signal sizes a hundredfold above the projection 
noise limit. 

PACS numbers: 42.62.Eh; 32.80.-t; 42.50.Pq; 42.50.Dv 



I. INTRODUCTION 



Optical atomic clocks based on neutral atoms confined 
in state independent optical lattices have made dramatic 
progress recently 0, [2, H[ . The highest spectral resolu- 
tion has been achieved in such a system [J], resulting in 
the clock instability approaching 1 x 10 -15 at 1 s and 
an overall uncertainty reaching 1 x 10~ 16 0, Q. The 
signal to noise ratio (SNR) for 10 4 atoms in these latest 
experiments is within a factor of two of the quantum pro- 
jection noise limit [6[. As the lattice clock performance 
approaches this limit squeezing of the collective atomic 
pseudo-spin to overcome quantum projection noise @, Q 
will lead to dramatic further advances in the clock per- 
formance because the number of atoms involved is large. 
Furthermore, we believe that the precision, control, and 
isolation from the environment achieved in these metro- 
logical lattice systems can be the basis of powerful probes 
to explore novel quantum dynamics, such as manifested 
in the collective interactions between an atomic ensemble 
and an optical cavity d M, HO, GJ, [H, [H, M, [U, HI • 

An atomic clock can be realized with the Ramsey tech- 
nique illustrated in Fig. [lja-c): atoms with two clock 
states \g) and |e) are driven with two 7r/2-pulses sepa- 
rated by a free evolution time T. The evolution of the 
atomic state during this clock sequence (Fig. D^b)) can 
be visualized on the Bloch sphere (Fig. [2(c)): The atomic 
pseudo spin points initially toward the south pole and is 
rotated around x by the first 7r/2-pulse to lie along y at 
position 1. From there the pseudo spin precesses along 
the equator to reach position 2 after T. This preces- 
sion is the "ticking" of the atomic clock: The total angle 
by which the pseudo-spin precesses measures the elapsed 
time. In order to measure this phase angle the position 
of the Bloch vector in the equatorial plane has to be 
translated into a position in the x-z-plane by the second 
7r/2-pulse to position 3 where it can be measured through 
the population difference between the levels which corre- 
sponds to the z-component of the pseudo-spin. 

The projection noise originates in the tip of the col- 
lective pseudo-spin not pointing in a sharp direction. 
Rather the position of the tip of the pseudo-spin is dis- 
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Figure 1: (Color online) (a) \g) and |e) are the two clock 
states for 87 Sr. (b) The Ramsey pulse sequence with each 
pulse having an area of 7r/2, and they are separated by a 
free evolution time T much longer than the duration of each 
pulse, (c) Illustration of the trajectory of the atomic pseudo 
spin on the Bloch sphere during the Ramsey pulse sequence. 
1: action of the first pulse, 2: free evolution, and 3: action 
of the second pulse. After the second pulse the population 
difference of the clock states is measured. 



tributed with a width of order V-/V for N independent 
atoms. The "hand" of the atomic clock is intrinsically 
fuzzy. As a consequence two phase angles closer to each 
other than 1/vN cannot be distinguished. This is il- 
lustrated in Fig. [2] where we show the sum of the prob- 
ability distributions of two Bloch vectors that we wish 
to distinguish. In Fig. [D(a) the two Bloch vectors have 
accumulated a relative phase of 1/VN at the end of the 
Ramsey pulse sequence at position 3 and hence the two 
Bloch vectors are not resolved. 

However, the projection noise can be reduced by 
preparing the atoms in an entangled state in which the 
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Figure 2: (Color online) Sum of probability distributions at 
the end of the clock sequence (position 3 in Fig. Q^c)) for two 
pseudo spin vectors that we wish to resolve in order to distin- 
guish two atom-field detuning values. The pseudo spins have 
accumulated a relative phase l/V~N during the free evolution 
stage of the Ramsey sequence. In (a) the initial state of the 
atoms was the conventional state of all atoms in \g) and the 
two final positions of the pseudo spin cannot be distinguished, 
giving rise to the shot-noise limit. In (b) on the other hand 
a spin squeezed state was used and the two positions can be 
resolved leading to a measurement precision better than the 
shot-noise limit. 



distribution of the pseudo spin is squeezed. The princi- 
ple of this idea is to prepare the atoms in a state such 
that the distribution of the pseudo spin at the end of the 
clock sequence is narrower in the z-direction, as shown 
in Fig. (Hb). This way the Bloch vectors of two states 
that have accumulated less than 1 / y/N of phase differ- 
ence can be distinguished. The higher phase resolution 
translates into an increased clock precision. By means 
of spin squeezing the projection noise can be reduced in 
principle to the Heisenberg limit which scales with the 
number of atoms as 1/N. The Heisenberg limit scaling 
has been demonstrated in experiments with few entan- 
gled ions [HI, [T^j . For neutral atoms, their large sample 
number permits a huge reduction of the projection noise 

tumi. 

Spin-squeezed states can be created by means of atom 
interactions [13] , interaction with a cavity field [HI, [24[ 
or through the back-action of quantum non-demolition 
(QND) measurements [2l|, [25|, (26|. Since we are consid- 
ering the latter approach in this paper we briefly review 
this type of measurement. For defmiteness we discuss 
the case of measurements of the z-component J z of the 
collective atomic pseudo spin. This is no restriction since 
squeezed states along any direction can be obtained from 
a state squeezed along the z-direction through rotations. 

According to the principles of quantum mechanics we 
are guaranteed to get the same result in a second mea- 
surement of any observable as in a first measurement of 
that observable provided that the observable commutes 
with the Hamiltonian. This can be considered an extreme 
case of squeezing: The probability distribution of the ob- 
servable has collapsed to a delta-function due to the first 
measurement. This type of measurement is sometimes 
called a projective measurement. The measurements that 



are of interest to us in this article are weaker than projec- 
tive measurements in the sense that several atomic states 
corresponding to different eigenvalues of the observable 
are consistent with the measurement outcome. The key 
point is that, conditioned on the measurement record, 
some states become more likely than others. This has 
qualitatively similar consequences for subsequent mea- 
surements of the observable as in the projective measure- 
ment case. As we keep measuring the observable we are 
more and more certain about the measurement outcome 
because only a small range of results will be consistent 
with the measurement record obtained up to that time. 
It is essential that these measurements be performed on 
the same system and not on different copies. Therefore 
the measurement has to be non-destructive. This type of 
measurement in which information about the state of the 
system is extracted gently in such a way that the state of 
the system persists after the measurement is complete is 
called a non-demolition measurement. As shown above 
the resulting state is squeezed. 

In our case the observable J z is a collective variable. 
Therefore J z cannot be measured through a measurement 
of each atom's spin, if we wish to achieve spin squeezing. 
That would project the atoms on a product state in which 
each atom's spin is pointing either up or down. Such a 
state has a reduced total angular momentum (J 2 ) and 
the reduction of the length of the angular momentum 
outweighs any possible benefits from squeezing. A mea- 
surement scheme for spin squeezing must therefore ensure 
that the total angular momentum J 2 is conserved. Math- 
ematically speaking this means that the Hamiltonian and 
the interactions describing the measurement must com- 
mute with J 2 . Physically it means that inhomogeneous 
broadening has to be negligible during the time scales of 
interest and the measurement must not be able to distin- 
guish between different particles. 

To summarize, a measurement protocol for the prepa- 
ration of spin-squeezed states should satisfy the following 
requirements. First, the probe must not lead to deco- 
herence by spontaneous emission or depolarization of the 
atomic sample by inhomogeneous effects (non-demolition 
and conservation of J 2 requirements). Second, the mea- 
surement must give the population difference with preci- 
sion exceeding the atomic projection noise. Finally, for 
clock applications it is important that the measurement 
does not introduce shifts of the clock transition. 

In this article we consider neutral atoms in an optical 
lattice clock. Conventionally the lattice is treated as an 
external potential in these systems, i.e. the back-action 
of the atoms on the light field is neglected. This approx- 
imation is motivated by the coupling between atoms and 
lattice photons being very weak. A large number of pho- 
tons is necessary to provide a deep enough lattice for the 
atoms while the many orders of magnitude smaller num- 
ber of atoms has only a microscopic effect on the light 
field. However, the minute changes that the lattice fields 
experience when they propagate through the atomic sam- 
ple contain information about the atomic state that is 
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normally lost. We show that the information that the 
atoms imprint on the lattice fields can be harnessed. In 
particular we propose to use the lattice field itself for the 
non-demolition measurement of the clock pseudo-spin to 
achieve spin squeezing. Such a scheme has several advan- 
tages over probing the atomic state with an additional 
interrogation field in addition to using a resource that is 
normally wasted. Importantly, the lattice does not intro- 
duce clock shifts as it operates at the magic frequency 
where the two clock states have an identical polarizabil- 
ity 0, 0, H, 0, H, El ■ Decoherence by spontaneous emis- 
sion is small since the lattice is far detuned from strong 
atomic transitions. The lattice laser also couples equally 
to atoms at different lattice sites due to the lattice peri- 
odicity, i.e. the probe is only sensitive to the total pseudo 
spin and the measurement is of the non-demolition type. 
According to our list of requirements above it remains 
to be shown that sufficient precision for spin squeezing 
can be achieved using this approach. That is the main 
subject of this article. 

In general one has to "level the playing field" between 
photons and atoms in order for the microscopic effects of 
the atoms on the light field to become detectable exper- 
imentally. Since it is difficult to increase the number of 
atoms this means that the number of photons has to be 
reduced. This can be achieved by putting the photons 
in a cavity. Figuratively speaking each photon passes 
through the atomic cloud many times and therefore a 
much smaller number of photons is sufficient to generate 
the optical lattice. Conversely each photon also accumu- 
lates the effect of interaction with the atoms over many 
round trips so that the signal is enhanced. In this article 
we consider the case of a bad cavity, in a sense that we 
will make precise below, for two reasons. First, this is the 
case that is immediately relevant for the next generation 
of atomic clocks. Second, we want to use the light field as 
a measurement device. This implies that the dynamics 
of the field should be as simple as possible so that the 
state of the atoms can be read off directly without hav- 
ing to understand the complicated physics of the meter. 
In a high-Q cavity where the atoms and light field inter- 
act with each other in the strong coupling regime their 
coupled dynamics would be so complicated that it would 
become hard if not impossible to infer the atomic state 
from measurements of the field. 

The principle of our idea for measuring the atoms' spin 
with the lattice fields is the following. The atoms are ini- 
tially prepared in state \g) . During the first Ramsey pulse 
the clock laser drives them into a 50/50 superposition. If 
this drive is in the Lamb-Dicke and resolved side-band 
regime the atoms will remain at rest. The recoil momen- 
tum associated with each transition has to be taken up 
by the lattice fields. If the lattice is generated in a ring 
cavity this is achieved by transferring photons from one 
mode to the counter propagating mode. By measuring 
the intensity redistribution one can determine how much 
momentum has been exchanged between the two modes 
or, equivalently, how many transitions have happened. 



This constitutes a measurement of J z since we started 
from a state with a known J z quantum number. 

The rest of this article is organized as follows. In sec- 
tion HQ we introduce the model and develop the theoret- 
ical framework that we will use. Section IIIII discusses 
the measurement scheme outlined in the previous para- 
graph. As we will see, the signal to noise ratio (SNR) of 
this measurement is insufficient to achieve spin squeez- 
ing for currently realizable lattice clocks. Section IIVI is 
dedicated to a superior measurement scheme based on 
detecting motional sidebands of the atoms that promises 
to lead to a strong enough signal that can yield significant 
spin squeezing. We draw conclusions in section fVl 

II. MODEL 

We consider N two- level atoms with ground state \g), 
excited state |e), and transition frequency uj a , trapped 
in a one-dimensional optical lattice generated by the two 
counter propagating running wave modes of a ring cav- 
ity with frequency ujl (Fig. [3]). The projections of 
the running waves' wave vectors along the z-axis are 
fei = cosQ lol/c and — where 9 is the angle at which 
clock laser and the lattice beams cross. The transverse 
profile of the modes is approximately Gaussian and we 
neglect the dependence of its radius wq on z. The cavity 
length is L. The atomic transition is probed by a highly 
stabilized clock laser of frequency uj c , which is linearly 
polarized in the same direction as the lattice. 

A. Effective Hamiltonian 

In our calculations, we neglect inhomogeneous effects. 
This is justified if the duration of the clock pulse sequence 
is smaller than the Ti time associated with the inhomo- 
geneities. The main source of inhomogeneous broaden- 
ing in experiments is due to the radial distribution of 
atoms in the ID optical lattice which is typically only 
weakly confining in the transverse direction. The inho- 
mogeneities stemming from this can be eliminated by 
confining the atoms in a 3D lattice. We assume that 
the wavefunction of the atoms can be factorized into a 
part for spin and z on the one hand and all other degrees 
of freedom on the other hand. Since we will be only in- 
terested in spin and z, considering just that part of the 
wavefunction is sufficient, leading to an effectively ID 
theory. The initial state, all atoms in \g) and vibrational 
ground state [28j . 

m = 0)) = \g,n = 0; g,n = 0;...; g,n = 0), (1) 

V v ' 

N times 

is completely symmetric under particle exchange with re- 
spect to the pseudo-spin and z which are the degrees of 
freedom relevant to us. The system will remain in the to- 
tally symmetric subspace since inhomogeneities are neg- 
ligible. Hence we can use a totally symmetrized basis 
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and going over to the usual second quantized formalism 
we can effectively treat the atoms as bosons. 



Figure 3: (Color online) Schematic of the lattice in the ring 
resonator with atoms trapped at the intensity anti-nodes. The 
clock laser propagates along the z axis and the projection of 
the wave numbers of the counter propagating lattice modes 
along the z-axis is fci and —kh- The lattice light leaking 
out of the resonator contains information about the atomic 
state. The clock laser wave vector is commensurate with the 
reciprocal lattice period under this particular lattice geome- 
try. Details of coupling into and out of the ring cavity as well 
as the detector arrangement are given in Fig. [4] 
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The system can be described with the Hamiltonian 



H = 



f dz i>l{z)^ a {z) + hg ]T /dz(e- tt **ak+e***al fc J(e^^ 

er— g,e z a—g,e 

+ ? ( e ~ M I dz e ^>e + H.c) + J2 i~^ala p + h V (al + a p )) . (2) 



ip e and ip g are bosonic field operators describing the 
atoms in the excited and ground states. The first term 
in the Hamiltonian is the atomic kinetic energy, with M 
being the atomic mass. The second term describes the 



lattice potential, go = 



is the coupling constant 



between atoms and the lattice field, a is the common po- 
larizability of \g) and |e). a±k L are bosonic field operators 
for the running wave modes of the ring resonator. The 
third term describes the drive of the atomic transition 
by the clock laser with a Rabi frequency O and detuning 
5 = uj c — u) a . The last term represents the detuning A of 
the lattice laser from the resonator mode and the pump 
of the lattice field with amplitude ~q, which is assumed to 
be real. The losses of both modes with decay constant k 
can be described by means of a standard Born-Markov 



Liouville operator 

p=±ki 



(al& p p ■ 



pal& p 



2a v pa\), (3) 



where p is the system density operator. We neglect spon- 
taneous emission from |e) which is justified because the 
lifetime of that level is many times longer than the clock 
sequence considered here. 

The periodic arrangement of atoms in the lattice 
strongly couples the two counter propagating modes with 
a coupling frequency of order goN. The symmetric co- 
sine and antisymmetric sine modes described by the field 
operators 



V2 



(a kL ±a-fcj 



(4) 
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on the other hand are uncoupled for a perfect atomic 
lattice and weakly coupled if the density distribution is 
slightly perturbed. Since it is exactly these small devi- 
ations of the atomic density distribution from the per- 
fect lattice caused by the clock laser that we wish to 
detect, we transform the lattice field to the symmetric 
and antisymmetric modes b + and 6_ . Because the cavity 
pumping amplitude rj is real only the symmetric mode 
is pumped. In steady state the symmetric mode 6+ has 
therefore a large mean amplitude /3+ which gives rise to 
the lattice. The steady state amplitude can be found 
from the equations of motion 



.d^) = _ {A + i K ){i+)+go f dz J2{TPl( Z )M Z )) X 

(2 cos 2 k L z (S+) + i sin 2k L z (£_)) + V2r], (5) 

where we have used that the expectation values like 
(& + ^>£Vv) etc. factorize in steady state to a good ap- 
proximation. We find /3 + by setting d(b + )/dt = 0, 



A + if -2 9o C(0)' 



(6) 



where we have used = in steady state and we have 
introduced 



C(t) = 



dz^cos 2 k L z(ipl{z)i}) a {z)). 



(7) 



We eliminate /?+ with the canonical transformation 6+ — > 
d + + /?+ , fe_ — > <i_ . /3+ gives rise to an optical lattice of 
depth 2h\[3 + \ 2 go. Assuming a deep lattice, we can neglect 
tunneling between different lattice sites and approximate 
the lattice with a harmonic potential at each site. 

The atom-lattice coupling is trivially lattice periodic. 
The clock-laser-atom interaction can also be made lat- 
tice periodic by crossing the lattice modes at an angle 



such that k c = 2/cl- Atoms on different sites are then 
indistinguishable and by transforming the coordinates of 
a particle trapped on site j according to z — > z — jft/kL 
we can treat them as if they were all trapped in a single 
harmonic well. This interpretation makes precise what is 
meant by the initial state Eq. ([1]). The indistinguisha- 
bility of atoms at different sites is also important from a 
practical point of view because it allows one to squeeze 
the pseudo-spin for all atoms, not just the atoms at each 
lattice site. 

We expand the atomic field operators in an energy ba- 
sis tfn(z) of the single harmonic oscillator representing 
the lattice, 



]fn{z)c a ,, 



(8) 



with oscillator frequency uj osc = y/4h\go\ \(3+\ 2 k 2 L /M and 
oscillator length a osc = yjhj (Muj osc ). We assume that 
the atoms are deep in the Lamb-Dicke regime which 
for atoms in the first few vibrational states means 

< 1. 

B. Mean field equations 



In the rest of this article, we study this system in the 
mean field approximation. The mean field equations of 
motion for the field amplitudes (d±) and the atomic am- 
plitudes {c a ,n) are found from the Liouville equation for 
the density matrix 



dp(t) 
dt 



[H,p]+C[p] 



(9) 



These equations close if we factorize the correlations be- 
tween atoms and light field according to (d±c a . n ) — > 
(4)<C ffl «) 1 (44) (4>(4), etc. We find 



.d(c e ,n) 
I : 

dt 



osc \^e,n/ 



n\e^ z \n>)(c g , n ,) 



(10) 



■d{Cg,n) 

dt 



d(d- 

dt 
d(d + 
dt 



+ihg ((/?; + (4»(4) - 08+ + (4»(4)) J2( n \ sin 2k L z\n'){c e , n ,) 



= nw OS c(cg,„) + — — 2j(n|e"' ifcc2: |n')(c e:n 

n' 

+thg ((/3; + (4»(4) - 08+ + (4))(4)) J2( n \ sin Mi.»\n')(& g , n >), 

n' 

= -(A - 2g S(t) + «c/2)(d_) - ig {{3+ + {d+))S 2 (t), 

= -(A-2g C(t)+iK/2)(L) + ig (d-)S 2 (t) + g {C 2 {t) - C 2 {Q))P+. 



(11) 



(12) 
(13) 
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We have introduced 



S(t) = ]T f dz sin 2 k L z(i>t(z)Mz)), (14) 
Sn(t) =J2[ dz sin 2fc L z(V-t (z)^(z)), (15) 



and 



C 2 (t) = J2f dz cos 2fcL2(V^(^)^(z)). 



(16) 



C. Approximations 

In the Lamb-Dicke and resolved sideband regime the 
density distribution of the atoms does not change much 
during the evolution. In particular we have 



C 2 (t)-C 2 (0) =0{{k L a 



) 2 )- 



(17) 



Neglecting this small term in Eq. (fT3|) is an excellent 
approximation for the cases we are interested in. 

With this approximation it is clear that the modes 
d- and d + are no longer pumped. Light is scattered 
into these modes exclusively through interaction with the 
atoms and we find the scalings 



9aS 2 {t) 



(d_), (d_) 



1% 



(18) 



As discussed in the introduction we consider the bad cav- 
ity limit. The previous equation motivates the appropri- 
ate condition for the bad cavity limit, 



9oS 2 {t) 



< 1. 



(19) 



The resulting hierarchy of the amplitudes (d+) <C 
(d-) <C /3+ allows us to make further approximations: 
We neglect the symmetric mode altogether, (d+) = 0, 
and we neglect the back-action of the field on the atoms 
in Eqs. (jlOllip which is of order g /3 + (d-) compared to 
the lattice potential of order <7o/3+ contained in lu osc . 

We end up with a theory in which the atoms move in 
the steady state lattice potential and are driven by the 
clock laser. Through S 2 (t) they are a source for the (d_) 
field. 



III. 



INTENSITY IMBALANCE 



As discussed in the introduction our goal is to find 
the population difference between electronic excited and 
ground states after a 7r/2-pulse by measuring the mo- 
mentum transfer between the +kt and — modes. A 
schematic of the detector arrangement for this measure- 
ment is shown in Fig. 0Ja). 



(a) 



a -k | resonator 



(b) 



| resonatoi 



S D- 



x/4 ^ 

^0 



Figure 4: In- and out-coupling into the ring cavity as well 
as detector arrangements for the two measurement schemes 
discussed in the text. In (a) two phase-locked lasers pump the 
two counter propagating modes of the ring cavity and the light 
leaking out of each mode is picked up with beam splitters. In 
(b) an incoming pump beam is sent through a Mach-Zehnder 
interferometer-like arrangement with a n/4 phase shift in one 
of the arms. This way the symmetric superposition 6+ of the 
two running wave modes is pumped. The d_ amplitude leaks 
out of the empty port of the Mach-Zehnder interferometer 
where its intensity is measured. 



To find the intensity imbalance between the modes we 
integrate the mean field equations of motion numerically. 
From the numerical solutions for the amplitudes inside 
the cavity we find the intensity difference at the cavity 
output 



51 = 



a-k L 



= 2nRe (3* + {L 



(20) 



Of course there will be noise in the measurement of the 
intensity imbalance. In order to find out whether this 
QND scheme is suitable for spin squeezing we have to 
determine whether it is possible to measure the number 
difference with a precision better than the atomic shot 
noise despite the noise in the intensity imbalance. 

We assume that the detection of the intensities of the 
two modes is (photon-) shot noise limited. The photon 

shot noise in each port is \l (a± k &±k L ) The result- 



ing SNR can be calculated as 



SNRi 



f™ dtK{al L a kL - aL kL a-kt) 
V 'itT dtn(al L a kL + a^a^J 



(21) 



Figure [5] shows the result as a function of the drive 
strength fl of the clock transition for realistic lattice 
clock parameters. The clock laser is resonant with the 
n = — * n = transition, i.e. 5 = 0. The lattice laser 
amplitude i] has been adjusted to give lu osc = 20u> rcci 
where ui rcc = huj j /(2Mc 2 ) is the recoil frequency at the 
lattice wavelength. We have assumed that the atoms are 
driven from ground to excited state with a 7r/2-pulse as 
they would during the first pulse of a Ramsey sequence. 
The detector outputs are recorded only during the pulse 
as indicated by formula Eq. (f2Tj) . 

In the resolved side band limit Q <C uj osc we can ob- 
tain an approximate analytical solution of the mean field 
equations of motion by assuming that the coherences be- 
tween different harmonic oscillator levels follow the drive 
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Figure 5: (Color online) SNR for detection of the intensity 
imbalance during a 7r/2-pulse for 10 6 87 Sr atoms in a ring 
cavity of length L = 1 cm, finesse F = 10 6 , and waist wq — 
30 fim. u osc = 20cu lac and 8 — (red solid line). The lattice 
laser is tuned to the magic wavelength Xl = 813nm and the 
polarizability at that wavelength is a — —5.37 x 10 _28 m 3 eo- 
The blue dashed line indicates the analytic result (Eq. J22}) 
in the adiabatic limit. 



adiabatically. From the adiabatic solutions for the atomic 
amplitudes we can then in turn find adiabatic solutions 
for the light field provided that the strength of the clock 
laser coupling is well within the band-width of the cav- 
ity, fl <C k. Once we have the amplitudes of the lattice 
modes we can calculate the signal to noise ratio using 
Eq. (EU). We find 



Nk c /k L 



16|pQ | Qujy, 



(22) 



As can be seen in Fig. the analytical solution agrees 
very well with the numerical results in the limit £1/lu osc <C 
1. 

The SNR falls off for small f2 because the signal 
strength is limited by a photon redistribution of order 
~ N/2 while the number of photons that contribute 
to the shot noise keeps increasing over longer intervals 
7r(2i7) _1 . At large fl the SNR falls off because more and 
more momentum is taken up by the atoms as they are 
driven harder and accordingly less intensity needs to be 
redistributed between the lattice modes to ensure con- 
servation of total momentum. The maximum SNR is 
obtained near Q rs lo osc [30] . Assuming that the expres- 
sion Eq. (|22p still holds, we estimate the maximum SNR 
as 



4 k r . 



SNR max = VN _j_ ^ /Wrec 
l-K kh V Wosc 



Ng 



(23) 



This simple scaling law is one of the key results of this pa- 
per. In order for the measurement to lead to spin squeez- 
ing, the SNR has to be sufficiently large such that the 
measurement uncertainty is smaller than the atomic pro- 
jection noise, i.e., SNR > N/ AN ~ y/~N. In other words 



the measurement has to be atomic shot noise limited. 
Since 4fc c (y / 7rfcL)~ 1 (w rcc /a; osc ) 1 / 2 is typically of order one 
and is hard to vary in experiments the collective coupling 
parameter NggK -1 is the all important parameter. Note 
that we can be in the strong collective coupling regime 
9} NgoK^ 1 ^> 1 while still in the bad cavity limit in the 
sense of Eq. (Tl9|) because 6>2 (t) <C N in the Lamb-Dickc 
and resolved side band regime. 

It can be seen from Eq. {23]) and Fig. [5] that a SNR 
greater than y/~N is hard to achieve with currently real- 
istic lattice clocks but it might not be completely out of 
reach in the future. The collective coupling parameter is 
small in these systems primarily because of the require- 
ment of operating the lattice at the magic wavelength. 
The magic wavelength is typically very far removed from 
atomic resonances. Therefore the atomic polarizability 
is very small, giving rise to a small single atom coupling 
constant go. The fundamental reason for the rather lim- 
ited SNR is the intrinsic photon shot noise of the lattice 
beams, i.e. the measurement is photon-shot noise lim- 
ited rather than atom-shot noise limited. The momen- 
tum transfer from the clock laser to the lattice laser that 
has to be measured is smaller than the large momentum 
uncertainties stemming from the photon shot noise in the 
lattice beams. In the conclusion section we discuss meth- 
ods that may improve the SNR of this scheme to a level 
where it becomes a viable means to obtain spin squeezing. 
However, we will first introduce a much more promising 
approach in the next section for measurement-induced 
spin squeezing. 



IV. SIDEBAND SPECTROSCOPY 

In this section we describe a superior detection method 
that is not affected by the large photon shot noise in the 
lattice beams that dominated the SNR in the detection 
method discussed in the previous section. The general 
idea is to design a measurement in which the signal can 
be separated from the lattice beams. The motional side- 
bands generated by the atoms oscillating in the lattice 
with frequencies ±o; osc are such a signal. In a hetero- 
dyne detector with bandwidth much smaller than uj osc 
the sidebands can be distinguished from the carrier at 
Wi, provided that the line width of the lattice is smaller 
than ct> sc which is achieved experimentally. 

In order for the detection of atomic vibration to con- 
stitute a measurement of the populations of atomic elec- 
tronic states these two quantities must be strongly cor- 
related. If for instance only atoms of one electronic state 
are oscillating while the other is at rest, detection of the 
sideband can measure the number of atoms in the oscil- 
lating state. Other possibilities would be to have atoms 
of both states oscillate with the same amplitude but n/2 
or 7r out of phase. If the atoms oscillate ir/2 out of phase 
the populations of the two states are proportional to the 
intensities in the two quadratures of the sidebands. If 
they oscillate it out of phase, the intensity of the side- 
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bands is directly proportional to the number difference 
of the two states. 

Such correlated states of atomic motion and internal 
level are created rather naturally in lattice clocks. To 
see this we consider the atomic dynamics during a ir/2- 
pulse starting from initial state Eq. ([1]). The limit of a 
strong pulse £1 3> lu sc is easy to understand intuitively. 
In this case the component of the atoms that undergoes 
a transition to the excited state receives N e Hk c of recoil 
momentum where N e is the final population of the ex- 
cited state, while the atoms that stay in the ground state 
remain at rest. After the 7r/2-pulse only the excited state 
component will oscillate at frequency o; osc . 

The atomic motion is more complex in the weak drive 
limit f2 -C lu osc . Because atoms exchange momentum 
with the lattice lasers while absorbing and reemitting 
clock photons it is clear that both states start oscillat- 
ing. The atoms' dynamics is illustrated in Fig. [SJ During 
the pulse the atoms in the two electronic states oscillate 
7r out of phase with each other but with equal envelope. 
The amplitude of the oscillations is suppressed compared 
to the strong pulse limit by a factor /c c a sc^/w sc- If 
the pulse duration is t-^/^osc = rn2iT with m an inte- 
ger, only the ground state oscillates after the pulse. For 
t?r /2^osc = (m + 1/2)2tt only the excited state oscillates 
and for t^/ 2 L0 O sc — {m ± 1/4)tt both states oscillate with 
equal amplitude. The oscillations after the pulse are 
undamped to a good approximation for k 3> lo OS c and 
A = 2g N regardless of the strength of f2 [H,[I3l- 

For the numerical examples in this section we consider 
a cavity with reduced finesse F = 10 4 (corresponding to 
k « 27T (3 MHz) for L = 1 cm) compared to the example 
in the previous section to ensure that n 3> uj osc . All other 
parameters are as before. 

We now turn to a quantitative discussion of the mea- 
surement of the atomic vibration using the lattice fields. 
As before we assume that we are in the bad cavity limit 
Eq. (fT5)) and we retain (<f_) and neglect (d+). We con- 
sider the situation outlined above where t^^Wosc — rain 
such that the excited state is at rest after the 7r/2 pulse. 
The amplitude is then fed by a source 



.9(1 



SNR 2 = 



2+T 



dt K \(d-)\ 2 



(25) 
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where T is the integration time and we start the mea- 
surement immediately after the tt/2- pulse. With a strong 
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Figure 6: (Color online) Momentum per atom in ground (a) 
and excited state (c) during a 7r/2-pulse with Q — 0.5w re c. 
Figures (b) and (d) show the corresponding phase space tra- 
jectories for atoms in ground and excited state, respectively. 
At the beginning of the pulse the atoms are at rest at the 
origin, p = and 2 = 0. During the pulse the atoms oscillate 
on ellipses of growing diameter. When the pulse terminates 
at t = 7r /2/fi the atoms in the ground state are at the turning 
point on the right. Therefore they continue to oscillate with 
large amplitude in the harmonic trap after the pulse. The 
excited state on the other hand is at the origin at the end 
of the pulse. It therefore remains essentially at rest after the 
pulse. 



local oscillator every d_ photon is detected and the SNR 
is equal to the total number of photons scattered from 
the symmetric mode into the aL mode. The SNR can be 
readily evaluated using the numerical solutions for the 
field amplitude. The result is shown in Fig. [7] for an in- 
tegration time of 1 s as a function of Q. As in the previous 
section the weak drive limit can be studied analytically 
by assuming that the coherences between the different 
atomic vibrational levels follow the drive adiabatically. 
The adiabatic result, 



SNRf ■ = 



{k c /k L ) 2 N 2 Q 2 K T 
32(A2 + («/2)2)|/3 + | 



(26) 



(n|sin 2k L z\n'){c\ l(J ){c n >^) w goN g z g sin(u; osc t-0 9 ) 

(24) 

where N g is the population of the ground state and z g 
and (f> g are the amplitude and phase of the oscillations of 
the ground state atoms. 

In order to evaluate the suitability of this scheme for 
spin squeezing we need to calculate the SNR for mea- 
suring the intensity of the sideband [31| . For het- 
erodyne detection with a strong local oscillator the SNR 
is 



is also shown in Fig. [7] and agrees well with the numer- 
ical result. The amplitude of the atoms' oscillations de- 
creases as the strength of the drive is reduced according 
to the aforementioned suppression by f2/ui osc , resulting 
in a weaker generated signal and accordingly a smaller 
SNR. 

In the strong drive limit the atoms eventually take up 
all the recoil momentum (N/2)hk c and the SNR saturates 
at 



SNR™ ax 



{k c k L a 2 osc ) 2 N 2 \P + \ 2 g 2 nT 
2(A 2 + (k/2) 2 ) 



(27) 



Figure [7] shows that for fl ~ ui rcc , a SNR >10 5 can be 
achieved. For 10 6 atoms, this corresponds to a measure- 
ment uncertainty a hundredfold smaller than the projec- 
tion noise, indicating that spin squeezing becomes pos- 
sible with this measurement scheme. Note that even at 
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Figure 7: (Color online) SNR for 1 s of detection at the vibra- 
tional sideband for 87 Sr. The asymptotic solutions for weak 
and strong drives, Eqs. (|26p and (|27[) . are shown as the dashed 
and dotted lines, respectively. Parameters are as in Fig. [5]ex- 
cept that the cavity finesse is F = 10 4 . 

fl ~ uJ rC c, the population of the first excited vibrational 
state is suppressed by a factor of (k c a OS c(Q / Uosc)) 2 ~ 
1CP 4 relative to the vibrational ground state, thus main- 
taining the validity of the Lamb-Dicke regime. 

V. CONCLUSION 

We have studied the feasibility of measuring the pop- 
ulations of clock states using the lattice field of an opti- 
cal lattice clock, with the goal of squeezing the atomic 
pseudo-spin. We have demonstrated that a measure- 
ment of population transfer based on detection of mo- 
mentum redistribution between the lattice beams is pos- 
sible but, with current experimental technology, not sen- 
sitive enough to enable spin squeezing. 

The reason for this failure is that the collective atom 
lattice coupling Ngo /n is too small for a lattice operating 
at the magic wavelength. A possible solution might be to 
operate the lattice at other magic frequencies where go is 
larger. Engineering the magic wavelength with external 
electric, magnetic or laser fields such that it is closer to 
an atomic resonance might be another option. Stronger 
coupling can also be achieved with a smaller volume cav- 
ity. 

A second measurement scheme based on detecting the 
motional sidebands of atoms after a clock pulse is promis- 
ing. The SNR that we calculated for realistic parameters 
appropriate for a 87 Sr lattice clock suggests that spin 



squeezing could be feasible with this method. 

Several questions will have to be addressed in the 
future. First it will be important to go beyond the 
harmonic approximation in which the optical lattice is 
treated as a single harmonic trap. Among other things 
this will allow us to elucidate the role of tunneling and the 
band structure during and after the clock pulses. Second 
the system needs to be studied in the good cavity limit. 
The good cavity limit is particularly interesting since our 
results indicate that the two detection schemes presented 
here work better in cavities with higher Q-factor. 

Already in the bad cavity, but especially in the good 
cavity case, it will be necessary to go beyond the mean 
field approximation. That would allow one to study in 
detail how the noise in one of the pseudo spin components 
is reduced at the expense of the other components. Such 
a beyond meanfield treatment is indispensable for finding 
the ultimate resolution limits of the spin measurements 
presented here. Possible routes are cumulant expansion 
of the atomic and lattice field operators, Langevin equa- 
tions and Monte-Carlo wavefunction methods. 

A question of fundamental interest arising in the 
context of the second measurement scheme is that in 
that case we are trying to measure states that are nei- 
ther eigenstates of the Hamiltonian nor of the operator 
s'm2kLZ with which the light field couples to the atoms. 
Instead we are trying to differentiate between states that 
differ in their dynamics, i.e. one of them is oscillating 
while the other is at rest. A measurement can there- 
fore not be done instantaneously since this would project 
the system on eigenstates of sin 2UlZ. Rather one has to 
carefully erase all information about the time at which a 
photon has been scattered from the symmetric mode into 
the antisymmetric mode by having the photons circulate 
in the cavity for a sufficiently long time before they leak 
out. 

We also plan to investigate the prospects of quantum 
feedback control [HI, H3] that should allow one to not 
only prepare the many-particle state probabilistically in 
squeezed states with a certain J z projection, but also to 
detcrministically drive the system to a target J z projec- 
tion. 
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